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A method is established for measuring low energy y-rays dose by using CMOS sensors without any X-/y-ray 
converters. Gamma-ray source of **'Am and "Eu are used to test the system. Based on gray value, an analysis 
method is proposed to obtain the y-ray dose. Cumulative dose is determined by correlating the gray value to the 
dose readings of standard dosimeters. The relationship between gray value and the cumulative dose of y-rays 
are trained by using back propagation neural network with BFGS algorithm. After comparison, it shows that 
BFGS algorithm trainings are suitable for different y-ray sources under higher error condition. These indicate 
the feasibility of measuring low energy y-ray dose by using common CMOS image sensors. 
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I. INTRODUCTION 


CMOS (complementary metal oxide semi-conductor) im- 
age sensor is advantageous in its low power consumption, low 
price and good radiation resistance [1-3]. In recent years, 
CMOS image sensors have been widely used in detecting X- 
and y-rays. Using a scintillating material as converters for 
X-ray to visible light photons, the CMOS image sensors are 
used in X-ray imaging detectors for medical imaging applica- 
tions [4]. However, the use of converters narrows the dynamic 
range [5]. The CMOS sensors are also used in y-ray detection 
and classification [6]. In this paper, we explore the possibil- 
ity of measuring low energy y-ray dose using CMOS sensor 
without any converters. Results of dosimetry measurements 
using 7+'Am and "Eu y-ray sources are presented and dis- 
cussed. The article is organized as follows. Fundamentals of 
low energy y-rays detection by using CMOS are described 
in Sec. II. In Sec. II, the gray value analysis method is in- 
troduced and fit the correlation between the energy of y-rays 
and the cumulative gray value, and a network model based 
on back propagation (BP) neural network and the Broyden- 
Fletcher-Goldfarb-Shanno (BFGS) algorithm is adopted to 
reduce the detection error and the convergence speed. A sum- 
mary is presented in Section IV. 


Il. METHODS 


The key parts used in this experiment are the UPA1021 
CMOS digital image sensors produced by Philips. The 
CMOS image sensor is composed of an image array of 
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Fig. 1. Schematic diagram of the system. 
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307200 (640 x 480) pixels, plus an on-pixel amplifier, tim- 
ing and control circuits, dual 8 bit A/D and video port, and 
signal processing unit [7]. The pixel size of the sensor is 
2.2um x 2.2 um, and optical size is 1/10 inches. The sen- 
sors are kept in black environment to reduce the interference 
by visible light in measurements. Besides, the CMOS image 
sensor is covered by an aluminum foil to prevent the « parti- 
cles. Fig. 1 is a schematic diagram of the system [8, 9]. 

The CMOS pixel unit is composed of a photodiode, an 
electrical signal conversion unit, a signal transmission transis- 
tors and a signal amplifier [10]. X- or y-ray photons produce 
photoelectric effect inside the photodiode. The photocurrent 
is proportional to the luminance component of CMOS. 

The digital signal processor (DSP) produces the driving 
signal for CMOS, and controls working condition of the 
CMOS. Also, it transforms the received signals into 2D im- 
ages, calculates the cumulative dose of y-rays, and sends 
the image and cumulative doses to the display screen. The 
SDRAM improves the image processing speed. 


A. Gray value analysis 


Gray value analysis is applied to analyze the results in this 
work. Gray value is brightness of an image or degree of color 
shades. Generally, gray value of pixel is acquired before com- 
plex image processing. Gray value analysis is widely used 
in satellite images, aerial photographs, geophysical observa- 
tions, etc. [11]. If the color values of a certain pixel in an 
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original image are (Red, Green, Blue), the gray value is 
calculated by [12] 


h = 0.299Red + 0.587Green + 0.114 Blue. (1) 


One y photon usually hits some adjacent pixels at the in- 
stant of imaging. Besides, crosstalk in CMOS image sensor 
always happens in adjacent pixels because of closeness of the 
pixels [13]. Then, cumulative gray value of CMOS pixels ir- 
radiated by the photons can be used to detect the y-rays. The 
gray value of one pixel is called as h(x, y), where (x, y) de- 
notes the coordinate of pixels around the brightest one in the 
light spot. The cumulative gray value can be calculated by: 


H =X h(z,y). (2) 


HI. BFGS ALGORITHM 


BP (back propagation) neural network is a multilayer feed 
forward network based on the error BP algorithm and is 
used to approach any continuous function in arbitrary pre- 
cision. Theoretically, a neural network with combinations 
of different training functions, nodes and layer numbers can 
fit any input and output nonlinear functions [14]. Many im- 
proved learning algorithm are proposed based on classical 
BP algorithm. In numerical optimization, BFGS is an itera- 
tive method to solve the unconstrained nonlinear optimization 
problems [15, 16]. As an improved the Newton’s method, it 
performs nicely even for non-smooth optimizations. In the 
MATLAB optimization toolbox, the BFGS algorithm is im- 
plemented by trainbfg [17]. 


IV. EXPERIMENTS 
A. The gray level response of “'Am 


The CMOS image sensors were irradiated by ~4!Am 
(1.26 x 10° Bq, 4.0534 x 10*y/s) source at Shanghai Insti- 
tute of Applied Physics, Chinese Academy of Science. In its 
œ decay nuclide, *41Am, emits 59.5 keV and 26keV y-rays. 
The CMOS chip was covered by an aluminum foil to prevent 
the « particles. 

First, the system was set in a radiation-free environment 
and the gray value response of *4!Am on CMOS is analyzed 
to test the background before the experiment. The results 
show the background gray value is 22. 

Next, the 2*!Am source was placed at 30 cm to the CMOS 
image sensor, and the y-rays were recorded at a rate of 30 fps. 
A typical image is shown in Fig. 2, where the pixel units 
not exposed by y-rays are kept black. The gray values of 
the lighted units were calculated by MATLAB, and plotted in 
Fig. 2. The light spot is composed of units with different gray 
values, some of which are lower than the background due to 
crosstalk in adjacent pixels. 
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Fig. 2. Typical image of a *'Am y photon acquired by CMOS and 
gray values of the lighted spots. 


B. Dose measurement 


The y photon response on the CMOS image sensor can 
be denoted by cumulating the gray values of all pixels. A 
standard dosimeter was used to record the radiation dose. The 
gray values obtained by CMOS and the standard doses were 
used to deduce the correlation. The dose on the dosimeter 
were recorded every three minutes, while the video of y-rays 
on CMOS was recorded continuously. Each frame of images 
was analyzed because the cumulative doses are denoted by 
the cumulative gray value of each y photon. 

The largest gray value of each frame was compared with 
the background. If the largest gray value exceeds background, 
the frame would be abandoned. Otherwise, the pixels of 
largest gray value were selected and the gray values of pix- 
els around the largest one were accumulated. The procedures 
were repeated with other images and the cumulative dose of 
241 Am was denoted by the accumulated gray value of all the 
images, which is called as H”. 

It was found that the cumulative gray values (H’) changes 
with sizes of the pixel matrix. To determine a proper size 
of the matrix, 13 x 13, 11 x 11 and 9 x 9 matrices were 
compared. A two-tier BP neural network (5 nodes in the in- 
put layer, 1 node in the output layer with the expected error 
1078) was created. The accumulated gray values (H’) were 
used as inputs, and the dose values recorded on the standard 
dosimeter were used as outputs. Thirty-two sets of sampling 
data were used to train the network by the BFGS algorithm. 

The correlation between the accumulated gray values (H’) 
and the cumulative dose for the *4!Am source is plotted in 
Fig. 3(a). The cumulative dose and accumulated gray value 
are in good correlation for all the three matrices. The H” is 
larger with a larger matrix under the same cumulative dose. 

The best method with proper matrix is selected by global 
consistency comparison. Data not used for training are often 
used as the testing data to estimate the global consistency, 
which is characterized by the dose error. Six sets of testing 
data were taken as the inputs of the trained network and the 
corresponding cumulative dose was the actual output. The 
actual dose (a) was acquired by the network and the expected 
dose (e) was measured by the standard dosimeter. Relative 
error of the dose (A) is calculated by 


A = [(a — e)/e] x 100%. (3) 
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Fig. 3. The correlation for *“'Am between the H’ and cumulative 
doses (a) and relative error of the cumulative dose (b), using the BP 
networks with BFGS trained by the data of 13 x 13, 11 x 11 and 
9 x 9 matrices. 


The dose errors using the three matrices are shown in 
Fig. 3(b). The dose error of the 9 x 9 matrix is the highest 
of all matrices, and the error of 13 x 13 matrix is the lowest 
of all matrices. Therefore, 13 x 13 matrix pixels are of better 
global convergence, and are chosen in the following analysis 
within the precision allowed. 


C. Comparison between different sources 


152Eu (9.97 kBq) is an important radionuclide for energy 
and efficiency calibration of y spectrometers [18]. '**Eu has 
a complex decay scheme, 27.9% by B7 emission and 72.1% 
by electron capture. It emits over 140 y-rays ranging from 
122keV to 1408keV. For comparison, '°*Eu y-rays were 
measured, with aluminum foil to block off the B particles. 


Nucl. Sci. Tech. 25, 060401 (2014) 


Relative error of 
cumulative dose (%) 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 
Accumulated gray value (x107) 


Fig. 4. (Color online) The testing dose error of the networks trained 
by data of 13 x 13 matrix of the '’Eu and ™'Am y sources. 


The same method was used to obtain the cumulative dose of 
152Eu, The 13 x 13 matrix pixels were used as the statistical. 
Testing dose error of the networks trained by data from the 
'S2Bu and 74!Am source are plotted in Fig. 4. 

It can be seen from Fig. 4 that the correlation between H’ 
and the cumulative dose trained for the *4!Am source by the 
BP network is also suitable for !*2Eu source, and the 24!Am 
has a lower testing dose error. The BFGS algorithm is em- 
ployed to compensate for the energy loss and reduce the mea- 
surement error in this experiment. Besides, the number of 
sampling points and the corresponding H” increase with the 
dose. This can offset some errors of the energy loss. There- 
fore, the error decreases as the accumulated gray value in- 
creases. 

From Fig. 4, the testing dose error of 74'Am differs from 
'S2Bu, due to the different response in CMOS pixels. In this 
work, the cumulative dose is measured with the known source 


type. 


D. Radiation damage analysis 


The performance (especially dark current) of CMOS is al- 
ways characterized by the gray value of dark images [19]. 
Radiation damage is analyzed by comparing the gray value of 
images acquired before and after the irradiation. The images 
were collected every three minutes. The largest gray value 
of each image was picked from pre-test and post-test of the 
experiment, and the results were both 22. After two days irra- 
diation by *4!Am and !°*Eu, only one dead pixel is found on 
the CMOS image sensor. The gray value of the dead pixel is 
less than 40. Experiment results showed that the dead pixels 
did not influence the dose error. Therefore, that exposure to 
small doses of radiation over short periods of time does not 
damage the CMOS or its functionality. 


V. SUMMARY 


In summary, a method is proposed to obtain the cumula- 
tive dose of **!Am and '*7Eu using common CMOS sensor 
without any converters, while the type of y source will not be 
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distinguished. The gray value analysis of the recorded images 
is adopted as the image preprocessing. The BP neural net- 
work and the BFGS algorithm are used to train the network 
and fit the correlation between the cumulative gray value and 
cumulative dose. Experimental results show that the y dose 
is calibrated in this method. The error of !°*Eu is lower than 
241 Am in the experiments. The simulation results show that 
the cumulative dose of low energy y-rays (7+!Am and !°*Eu) 
detection is completed. The irradiation experiments show that 
the irradiation damage is very little and can be ignored. Based 
on the experimental results, it is proposed that the general 
CMOS image sensors can be used to measure the cumulative 
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